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Abstract. A new class of high temperature superconductors based on iron and 
arsenic was recently discovered [T], with superconducting transition temperature as 
high as 55 K [5] . Here we show, using microscopic theory, that the normal state of the 
iron pnictides at high temperatures is highly anomalous, displaying a very enhanced 
magnetic susceptibility and a linear temperature dependence of the resistivity. Below 
a coherence scale T*, the resistivity sharply drops and susceptibility crosses over to 
Pauli-like temperature dependence. Remarkably, the coherence-incoherence crossover 
temperature is a very strong function of the strength of the Hund's rule coupling Jnund- 
On the basis of the normal state properties, we estimate Juund to be ~ 0.35 eV. 
In the atomic limit, this value of Jnund leads to the critical ratio of the exchange 
constants J1/J2 ~ 2. While normal state incoherence is in common to all strongly 
correlated superconductors, the mechanism for emergence of the incoherent state in 
iron-oxypnictides, is unique due to its multiorbital electronic structure. 
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The unusually high superconducting critical temperatures in iron-oxypnictides 
together with unusual normal state properties, which do not fit within the standard 
framework of the Fermi liquid theory of solids, place the iron pnictides in the broad 
category of strongly correlated superconductors, such as k organics, cerium and 
plutonium based heavy fermions, and cuprate high temperature superconductors. In all 
these materials, superconductivity emerges in close proximity to an incoherent state with 
unconventional spin dynamics, that cannot be describe in terms of weakly interacting 
quasiparticles. Describing the normal state of these materials, is one of the grand 
challenges in condensed matter theory and has resulted in numerous controversies in 
the context of the cuprates. 

Iron pnictides show very high resistivity and very large uniform susceptibility in the 
normal state [T]. Superconductivity emerges from a state of matter with highly enhanced 
Pauli like type of susceptibility [3]. As a function of F~ doping, the susceptibility is 
peaked around 5% doping reaching a value 25 times bigger than Pauli susceptibility 
given by LDA [Hill]- In the parent compound, the resistivity exhibits a peak at 150 K [1] 
followed by a sharp drop, which is due to a structural transition from tetragonal (space 
group P4/nmm) to orthorhombic structure (space group Cmma) E] followed by a 
spin density wave transition [51 [6j at lower temperature. Only a 5% doping completely 
suppress the specific heat anomaly jH], while resistivity still shows a very steep drop 
below a characteristic temperature T* P [9] . While there are suggestions that the drop 
of resistivity is due to opening of a pseudogap and proximity to a quantum critical 
point in samarium compound [9j, other measurements in lanthanum compound seem to 
suggest less exotic and more Fermi liquid-like resistivity at 10% doping level [T0~t[TT]. Here 
we will show that the steep drop of resistivity in doped compounds can be understood on 
the basis of incoherence-coherence crossover, which is naturally accompanied by gradual 
screening of the magnetic moments, leading to enhanced magnetic susceptibility, and 
crossover from Currie- Weiss to Pauli type susceptibility. 

Theoretically, electron-phonon coupling is not sufficient to explain superconductiv- 
ity in LaOi_ x F x FeAs [TSJ [32] and spin and orbital fluctuations need to be considered. 
The weak coupling approach based on RPA approximation has been carried out by 
several groups for a simplified two band model containing two orbitals, d xz and d yz 
[TU HHJ [16]. The nature of the superconducting state and the pairing symmetry is 
apparently very sensitively to the parameters of the model, and therefore it is very 
desirable to reliably determine parameters of a microscopic Hamiltonian, which is the 
purpose of this article. 

The electronic structure of the LaOi-^F^FeAs compound was studied by the density 
functional theory (DFT) |3] and the dynamical mean field theory (DMFT) [12]. DFT 
calculations show that this material has a well separated set of bands with mostly d 
character in the interval between —2 eV and 2 eV around the Fermi level [12]. The parent 
compound has five Fermi surface sheets [I], two electron cylinders, centered around the 
zone edge (M-A line), and three hole pockets around the zone center (r~Z line). Upon 
doping with electrons, the hole-like pocket quickly disappear. The nesting wave vector 
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between the hole and electron pocket promotes the spin density wave instability which 
was actually observed in experiment [6] . Notice however that the Bragg peak in neutron 
experiments remains commensurate and does not change with doping, pointing to the 
inadequacy of the weak coupling spin density wave scenario. 
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Figure 1. a) Sketch of the orbital levels in the downfolded LDA Hamiltonian for 



LaOFeAs, keeping five d-orbitals (see appendix Appendix A I. The levels are obtained 
by h a 



LDA 

k,QCK ' 



where a is the orbital index. The Hund's coupling in this 



compound becomes relevant when its strength is comparable to the total splitting 
J ~ 350 meV. The coordinate system is chosen such that x and y axis point from Fe 
atom towards its the nearest neighbor Fe atoms, b) Probability for each iron 3d atomic 
state in DMFT calculation for LaOFeAs at T = 116 K and J Hun d = 0.4 eV. There are 
1024 atomic states in Fe-3d shell. Within sector with constant occupancy states are 
sorted by increasing atomic energy). 



Fig. [T] sketches the crystal field levels, as obtained by the Local Density 
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Approximation (LDA), when the localized Wannier orbitals are constructed for Fe- 
3d states. Notice that the tetragonal splittings are comparable to the t2g-eg splittings, 
leaving only two degenerate orbitals, i.e., xz and yz. Notice however, that the bandwidth 
of these orbitals, as shown in Fig. 3 of Ref. [12j, is of the order of 3eV, which is one order 
of magnitude larger than the crystal field splittings. 

A naive atomic picture suggests a large magnetic moment on the Fe atom (S = 2), 
and consequently almost classical moments, which can not condense by Copper pairing, 
but would rather order magnetically, just like they do in pure iron. Due to degeneracy 
of xz and yz orbital, even an infinitesimal Hund's rule coupling would therefore lead to 
a spin state of at least 5 = 1. 

Dynamical mean field studies [12], showed that the correlations in the parent 
compound are strong enough to make this material a bad metal with a large electronic 
scattering rate at, and below room temperature, as seen in the momentum resolved 
photoemission spectra, and the absence of a Drude peak in the optical conductivity. 
However, it was argued in Ref. [12], that the correlations are not strong enough to open 
the Mott gap in the parent compound, in contrast to the hole doped cuprates, which 
are known to be doped Mott insulators. The DMFT study determined LaOFeAs, to be 
a strongly correlated metal, which is not well described by either atomic physics, nor 
band theory. Furthermore, it was shown that all the five iron d orbitals have appreciable 
one particle density of states within 0.5 eV of the Fermi level, and therefore a realistic 
minimal model for these materials should contain all five d iron orbitals. 

In this article, we use the LDA+DMFT method p2] in combination with the 
continuous time Quantum Monte Carlo method [18], to investigate the transport and 
thermodynamic properties of the 10% doped LaOi_ x F x FeAs compound. More details 
of the method are given in Ref. [T2l [T7] . 

Our goal is to understand why a material with high orbital degeneracy, which 
naively would have expected to be an excellent metal [19], exhibits strong correlation 
effects, and what determines the crossover from the incoherent regime to the coherent 
Fermi liquid state. We show that the Hund's rule coupling Jnund dramatically reduces 
the coherence scale T* in these compounds, and promotes the highly incoherent metallic 
state. 

In physical terms, correlated quasiparticles are spin one half objects that have to 
be built from atomic states, which are dominantly S = 2 and S = 1, in the presence of 
Hund's rule coupling. As a result, the overlap of the quasiparticle and the bare electronic 
states is very small. 

To understand the nature of the metallic incoherent state, it is useful to study the 
probability that an electron in this compound is found in any of the atomic states of 
iron 3d orbital. This is plotted in the atomic histogram in Fig. [I] introduced in Ref. [20] . 
Even the most probable atomic states have probability of only a few percent, hence 
a naive atomic limit is qualitatively wrong for this compound. The atomic histogram, 
characterizes the nature of the ground state of the material. To investigate the coherence 
incoherence crossover, we need to probe the excitation spectra, through the transport 
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and the thermodynamic probes. 
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Figure 2. a) The local susceptibility versus temperature for few different Hund's 
coupling strengths, b) Resistivity versus temperature for the same Hund's coupling 
strengths. 



We first turn to the magnetic susceptibility which is the most direct probe of the 
spin response: 

k B Jo 

Here S z is the spin of the iron atom, is a magnetic moment of a free electron, Na 

is the Avogadro number, and f3 is the inverse temperature. In the absence of screening 
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(static S z ), susceptibility should follow the Curie- Weiss law \ — ^ A ^ s+V > , In the 
opposite itinerant limit, the spin susceptibility exhibits Pauli response \ — Hb-NaDq, 
where Do is the density of state at zero frequency. In Fig. [2^, we show the local spin 
susceptibility for a few values of Hund's coupling Jnund in 10% doped compound. 

In the absence of JHund, the susceptibility obeys the Pauli law, and the mass 
enhancement due to correlations is negligible. On the other hand, at large JHund = 
0.7 eV, the system remains in the local moment regime down to lowest temperature we 
reached (50 K). For some intermediate values of Hund's coupling Juund ~ 0.35 — 0.4 eV, 
the system evolves from an incoherent metal at high temperature, to the coherent Fermi 
liquid at low temperature. 

The metallicity of the compound always leads to screening of the electron magnetic 
moments, which eventually disappear at sufficiently low temperature and frequencies, 
but are present at high temperature and high energies. The scale at which the screening 
of the local moment takes place, however, is a very sensitive function of of the strength 
of the Hund's coupling, being very large for Juund = and unobservably small for 
JHund = 0.7 in the temperature range considered. This is one of the main messages of 
this article. 

To shed light into the coherence-incoherence crossover, we computed the resistivity 
due to electron-electron interactions, given by 

\ = Vohl dU (~<£) II Tr M w )Pk(^HMPkM] (2) 

where are velocities of electrons, pk is the spectral density of the electrons, and V is 
the volume of the unit cell. The resistivity in Fig. [2]d demonstrates the crossover even 
clearly. The resistivity is approximately linear with modest slope above the coherence 
crossover temperature, and it exhibits a steep drop around the coherence temperature. 
It is important to stress that the drop of resistivity here is not due to proximity to a 
quantum critical point |9J, or due to proximity to spin density wave j5], but because 
the electrons, being rather localized at high temperature, start to form a coherent 
quasiparticle bands with the underlying Fermi surface. 

In most compounds, the strength of the correlations is controlled by the ratio of 
the diagonal on-site Coulomb interaction U to the bandwidth, while Hund's coupling 
usually plays a sub-leading role. We find that in this compound, a crucial role is play 
by the Hund's coupling. This point was first noticed in the context of a two band 
Hubbard model calculation at half-filling [2T] , and it was shown that the coherence scale 
is exponentially suppressed when Hund's coupling is taken into account. Moreover, it 
was emphasized that the non-rot at ionally invariant Hund's coupling, usually assumed 
in combination with traditional quantum Monte Carlo methods [22], can lead to a 
qualitatively and quantitatively wrong results. Indeed, we have found that the absence 
of the rotational invariant Hund's coupling leads to a first order phase transition 
as a function of JHund, between a Fermi liquid and non-Fermi liquid state at zero 
temperature. On the other hand, the rotational invariance of the Coulomb interaction 
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recovers Fermi liquid at zero temperature, although the coherence scale can become very 
small, as seen in Fig. ^jp for Jnund = 0.7. 

It is interesting to note that the atomic Jnund in iron is approximately ~ 1.2 eV. 
In many materials with correlated / orbitals, the screening of Jnund in solid is well 
accounted for by 20% reduction from its atomic value. In LaOFeAs compound, only 30% 
of the atomic Hund's coupling gives the best agreement with experiment for resistivity 
and susceptibility. 

To elucidate the sensitivity of mass enhancement with varying Hund's coupling, 
we show low temperature specific heat in Fig. [3j For small JHund, we find specific 
heat value very close to the LDA prediction 7 ~ 6mJ/(molK 2 ). For Jnund = 0.35 eV, 
a value which gives resistivity in very favorable agreement with experiments, requires 
the enhancement of the specific heat for factor of 5 to 30 m J /(molK 2 ). This is large 
enhancement and should be easy to measurable in future experiments, and is a key 
prediction of the theory. This would be a strong test that a local theory, as the one 
described in this paper, correctly predicts the physical properties of this material, since 
would allow to differentiate it from paramagnon or itinerant spin density wave theories, 
where the enhancement of the static susceptibility is much larger than the enhancement 
of the specific heat. At this point, reliable data for the normal state specific heat of the 
doped compound is not available down to sufficiently low temperature to reliably extract 
the coefficient 7. It would be interesting to measure it in the slightly doped compound 
(5%), where the superconducting temperature is low enough, and susceptibility shows 
the strongest enhancement. 

The full computation of the momentum dependent spin response within 
LDA+DMFT is still a formidable task. However one can gain insights into the 
momentum dependence of \ by making a simplifying assumptions, that near the Mott 
transition xIq)" 1 = a + Jq, where J q is the Fourier transform of the exchange constants, 
and a is a local but frequency dependent object. 

We computed the first neighbor J\ and next-nearest neighbor J2 exchange constants 
as a function of Hund's coupling. For this purpose, we first constructed Wannier orbitals 
from iron 3d bands (see appendix Appendix A), we exactly diagonalized the iron atom, 
and perform the second order perturbation theory, with respect to iron-iron hopping. 
This perturbation is carried out around the atomic ground state of the S — 2 and 5=1 
sector. 

Fig |HJd shows the exchange constants as a function of the Hund's rule coupling. 
This provides complementary information to the exchange constants computed around 
the band limit, which were recently determined by LDA [2"H [231 [25] , and are typically 
larger in magnitude than those obtained here. The exchange constants are supportive 
of the actual ordering in the parent compound, determined recently by the neutron 
scattering experiment [7] , if J2 > Ji/2. For reasonable Hund's coupling, both S — 1 
and 5* = 2 ground states respect this inequality. If we select the S = 1 state for 
the ground state, the next-nearest J 2 is much bigger than J 1 . For the actual ground 
state of the atom (S = 2), and realistic Jnund ~ 0.35, the value of J2 is close to 
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Figure 3. a) Low temperature Fermi liquid specific heat coefficient 7 versus 
Hund's coupling strength. The effective mass is exponentially enhanced by the 
Hund's coupling. The coherence scale is consequentially exponentially suppressed 
by the Hund's coupling, b) Exchange constants, as obtained from the second order 
perturbation theory around the atomic limit, considering either S=2,N=6 or S=1,N=6 
for the ground state. Positive values for the exchange constants correspond to 
antifcrromagnetic couplings. The rest of the exchange constants are small, for example, 
J 3 ~ 1.8K and J z ~ 0.2 K. 



J 2 >= Ji/2. It is interesting to note that according to a well known result [28], this 
particular choice of exchange constants results in large magnetic fluctuations due to 
frustration. At temperature above the spin density wave transition, there can also be a 
second order Ising phase transition, at which the two sublattice magnetizations of the 
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two interpenetrating Neel sublattices, lock together. According to Ref. [28], the Ising 
transition would occur around T c ~ J 2 S 2 ~ 200K. 

The large value of the next nearest neighbor exchange constant demonstrates that 
even in the atomic limit, the system is highly frustrated, as pointed out in Ref. [2U [27J, 
which gives us more confidence that the single site DMFT description, employed in 
this article, is adequate, and the momentum dependent self-energy, as obtained by the 
cluster extensions of DMFT [T7], will not qualitatively change the physics. 

From the above provided evidence, LaOi_ x F x FeAs system fits nicely within 
the general rubric of strongly correlated superconductivity. The hallmark of this 
phenomena, is the appearance of a incoherent normal state, above T c . In the classic 
Bardeen-Cooper-Schrieffer theory of superconductivity, the effective Fermi energy is 
much larger than the pairing interaction, resulting in small critical temperatures. Strong 
correlations play the role of reducing the the kinetic energy, allowing the pairing 
interactions to lead to appreciable critical temperatures |29j. 

While these are common themes shared by all the high T c materials, the mechanisms 
by which different systems manage to renormalize down their kinetic energies, and the 
details of the competing states that appear as a result, are different in different classes 
of materials. 

Hole doped cuprate superconductors can be modeled as a single-band materials, 
with Coulomb correlation slightly above the Mott transition (in DMFT terminology 
the effective U, which arises from the charge transfer energy, is above U C 2), and indeed 
maximal T c is achieved when U is in the vicinity of critical U C 2- 

We have argued in our earlier publication [12], that the LaOFeAs is an intrinsically 
multiorbital system, which is slightly below U C 2- The mechanism which is responsible 
for the incoherence of the normal state and the renormalization of the kinetic energy is 
the Hund's rule Jnund, which we estimate here to be Juund = 0.35 — 0.4 eV. Competing 
states in this material are magnetic states, that have been predicted in Ref. [5] and since 
observed by neutron scattering experiments [7J EJ |3]. 

The incoherent normal states, is then eliminated by either a magnetic or a 
superconducting state, and the study of the competition between these mechanism in 
the multiorbital framework, deserves further attention. 



Appendix A. The tight binding Hamiltonian for the Fe-As planes 

The Fe-As planes in LaFeAsO can be described by a tight-binding Hamiltonian 



containing Fe and As atoms, with two Fe and two As per unit cell (see Fig. Al 
for drawing of the unit cell containing Fei, Fe2, Asi and AS2). We obtained the 
tight-binding parameterization by downfolding the Local Density Approximation first- 
principles results, which are expressed in a well localized Muffin Tin Orbital (MTO) base. 
The localized nature of the MTO base makes the projection of Kohn-Sham orbitals to 
a localized orbital basis set simple, and downfolding procedure more constrained and 
well defined. 
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E - onsite LaFeAsO LaFePO 

E[x 2 - y 2 ] -0.510 -0.654 

E[yz] -0.535 -0.656 

E[z 2 } -0.564 -0.635 

E[xy) -0.650 -0.692 

E[ Pz ] -2.144 -1.990 

E]p x ] -2.440 -2.358 



Table Al. On-site energies for both LaFeAsO and LaFePO in the model with both 
Fe-d and As-p orbitals. 



The few largest hoppings are sketched in Fig. AT and A2 All energies and hoppings 
are in units of eV. The on-site energies of all orbitals are tabulated in table [AT] and most 



of the hoppings are tabulated in tables A2 and A3 Note that orbitals are expressed 
in the coordinate system of the four atoms unit cell, shown in Fig. Al The complete 
set of tight binding hoppings is available to download |23J. For comparison we also 
include the hoppings for a closely related compound LaFePO, which shows considerably 
less correlation effects due to somewhat larger crystal field splittings, which renormalize 
down the effective Hund's coupling. 




Figure Al. left top) choice of the unit cell in the FeAs plane. Asi lays above the Fe 
plane, and As2 below it. 1-4) Sketch of some largest hoppings and their values. 
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Figure A2. Sketch of a few more relevant Fe-As hoppings. 



Further downfolding to the 5 band model of Fe is possible. The on-site energies and 
the few largest hoppings of the corresponding multiorbital Hubbard model are tabulated 
in table 3. The sketch of the most relevant hoppings is shown in Fig. X3 Note that the 
coordinate system here is chosen differently than above. Here we use the 45° rotated 
coordinate system, which corresponds to usual square lattice of irons. 
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hopping figure LaFeAsO LaFePO 
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Table A2. The few largest hopping matrix elements in the model of Fe-d and As-p 
orbitals. Each Fe atom is surrounded by two incquivalent As-atoms. The two different 
As orbitals are sketched in Fig. A1|A2 and are denoted by p and p'. 
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Table A3. Continued from table |A2[ a few more large hopping matrix elements of 
the Fe-d,As-p model. 
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